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The zero gap surface states of a 3D-topological insulator host Dirac fermions with spin locked 
to the momentum. The gap-less Dirac fermions exhibit electronic behaviour different from those 
predicted in conventional materials. While calculations based on a simple linear dispersion can 
account for observed experimental patterns, a more accurate match is obtained by including higher 
order k terms in the Hamiltonian. In this work, in presence of a time reversal symmetry breaking 
external magnetic field and higher order warping term, alteration to the topologically ordained Berry 
phase of (2n+ l)7r, momentum relaxation time, and the magneto-conductivity tensors is established. 


I. Introduction 

Topological insulators (TI) such as 612863 and 612 X 63 
ar 6 charact 6 riz 6 d by surfac 6 stat 6 S which possess distinct 
characteristics, for instance, a non-trivial 6 erry phase of 
(2n + l) 7 r, mass-less Dirac fermions, and a spin helical 
structure under ideal conditionsTime reversal sym- 
metry(TRS) breaking perturbations destroy the topolog¬ 
ical nature of the surface states and introduce additional 
contributions that have been experimentally established 
primarily through ARPES data4 Further, the linear dis¬ 
persion, at points away from the Dirac cone is altered by 
higher order fc-terms^ that in turn modify the topologi¬ 
cally ordained behaviour.— An external magnetic field 
which breaks TRS and higher order fc-terms that warp 
the band structure is considered in this work to 1 ) com¬ 
pute deviations to the Berry phase from the well-known 
( 2 n-|- l) 7 r and 2 ) evaluate the momentum relaxation time 
for electrons scattered on the surface of a 3D-topological 
insulator. 

The paper is outlined as follows: In Section II, an ex¬ 
pression for the Berry phase using dispersion and wave 
functions obtained from the warped and a gapped TI 
Hamiltonian is calculated which explicitly shows a non- 
topological component followed by a determination of the 
momentum relaxation time within a Boltzmann equation 
framework. The expressions derived in this section are 
numerically evaluated in Section III. A brief summary 
of results and their potential implications concludes the 
paper. 

II. Theory and Model 

The two-dimensional Dirac Hamiltonian is^di 

Hsurf.states — ^Vfi^CTxky fJykx) ( 1 ) 

where Vf denotes Fermi-velocity and ai;i = x,y are the 
Pauli matrices. The two-dimensional Dirac Hamiltonian 
is in principle sufficient to probe the surface states; how¬ 
ever, ARPES studies of the Fermi-svxiace at energies sig¬ 
nificantly far away from the Dirac point reveals a snow¬ 
flake like hexagram structureii^— that is markedly dif¬ 
ferent from a simple circular Fermi surface observed by 
application of Eq. [T] This departure from experiment 
is reconciled by noting that the simple two-dimensional 


Dirac Hamiltonian fails to account for the underlying 
crystal symmetries. The deformation of the Fermi sur¬ 
face can be theoretically reproduced if higher order k 
terms are incorporated in the Hamiltonian. Since the 
two-dimensional Dirac Hamiltonian must comply with 
the C 3 „ point-group and time reversal symmetry, the 
next higher order terms that must be added are cubic 
in k. The modified Hamiltoniani^ therefore must look 
like 

H{k) = eoik)+hvf{axky-aykx) + ^h^ {k^ + k^) (2) 

eo{k) introduces the particle-hole anisotropy and the cu¬ 
bic terms denote warping. Using Eq. [21 and ignoring 
particle-hole anisotropy, the surface state spectrum is 

e±(fc) = eo{k) ± hv'jk'^ + X^hfk^cos^{W) (3) 

where 9 = tan~^ 1 ~ ■ The spectrum contains the low- 

\kx J 

est order correction to the perfect helicity of the Dirac 
cone predicted in Eq. |TJ The cos'^{2>9) term possesses the 
symmetry of the C 3 ^ point group and the Hamiltonian is 
evidently time reversal symmetric. 

A. Berry phase of gapped Dirac fermions 
The Berry phase of the gapped surface spectrum of a 
3D topological insulator was derived in Ref. [H. A finite 
band gap can be induced either through the proximity 
effect of a ferromagnet or an s-wave superconductor— 
The final expression for Berry phase is given as 


7^ = TT 


I 1+ - 

\ \J ^lro + {kVfkf 


(4) 


where Ap^o is the band gap split on account of proximity 
effects and r; = ±1 denotes the helicity. The goal of this 
section is to examine the influence of warping on the final 
Berry phase. To begin, the wave functions (I'll (r;i?))) 
of the warped Hamiltonian inserted in the Schrddinger 
equation (r;i?)) = U„(i?)|'I' (r;i?)) can be writ¬ 

ten as 




1 
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f Xri{k)exp{—i9)\ 

V ) 


(5a) 
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where 


A„(fc) = 


2 ^ ^ _ ^pro _ 

\ \l ^'pro + 


(5b) 


the usual relations kx = kcos{0) and ky = ksin{9) are 
modified to k^ = k{9)cos{9) and ky = k{9)sin{9). The 
angular derivatives are therefore 


dkx dk n 1 ■ n 

—r-- = —cos9 — ksin9 
d9 d9 


(12a) 


^pro ~ ^pro "b 2 (^+ + ^-) 


(5c) 


The warping term effectively augments the band gap 
splitting originally introduced in the topological insu¬ 
lator. The complete Berry phase in fc-space such that 
R = k can be written as 


7 „(C') =*£(T(r;R)|Vfl|^'(r;i?))dR 


( 6 ) 


Expanding the Berry connection Ay = 
(d^ (r; i?) |V_r| 4'(r; i?)) in Eq. [6] using the wave 
functions given in Eq. [S^ yields 


= 1 [X*exp{i9) r]X*_y) (^\dk9) exp{ i9) 


pdkX- 


Simplifying the above expression, 

Mk) = '^{X;d,Xy-i\Xy\^d,9 + X%dkX-y) 
= l\Xy{k)\^dk9 


(7) 


( 8 ) 


The last expression has been condensed by noting that 
dkX-y evaluates exactly as dkXy with sign reversed, 
therefore taken together they are equal to zero. dkXy 
is worked out. 


dkXy(k) = dk 1 ± ,— 

\ \/^ 




pro 


{hvfk) 


Vdk\J Ap^o + {kvfkf 


(9) 


The final expression for Berry connection in matrix no¬ 
tation is therefore 


= 2 I 1 ± 


A„ 


1 f-K 


\J ^pro + {kVfk) 


2 / 


( 10 ) 


In deriving Eq. [TOl the angular derivatives dk^9{k) = 
k k 

—and dky9{k) = were used. Inserting Eg. fTOl in 
Eq. [6] and integrating over a closed path in fc-space gives 


and 


dky dk 

= —sin9 + kcos9 
d9 d9 


(12b) 


The Berry phase with warping using Eq. [TT] and substi¬ 
tuting for angular derivatives from Eq. [12] gives 


1 / I Apro \ -ksin9 


I 1 ± 


A, 


pro 


kcos9 


r 2 Tr 


— r, 


d9 \l± 


\J Apio + {hvfkf j 

A -I- Xh^k^cos{W) 

\J {A + Xh^k^ cos{i9yf + {hvfk)^ 

(13) 


An analytic evaluation of the integral is not possible and 
a numerical solution is presented in Section III. The in¬ 
tegral, as evident, consists of the topological contribu¬ 
tion of TT and the non-topological part shown in Eq. 1141 
The warping term effectively increases the proximity in¬ 
duced band gap and the non-topological component of 
the Berry phase. 


“^non-top — 


A -I- Xh^k^cos{W) 


\J {A + Xh^k^cos{Wyf + {hvfkY 


(14) 


The Berry phase induced fictitious magnetic fiel d^^i^° by 
computing y x and setting the warping term to zero 
is 


>/ic — T; 


1 h^v'iA 


y A-ipro 


+ (kVfky 


(15) 


The additional fictitious magnetic field generates a veloc- 

• • dk . ... 

ity given as Va = — x B jic- This velocity is transverse 

to the electric field and gives rise to the intrinsic Hall 

current. Writing — = —-^E, the corresponding Hall 
dt h 

current is 


= (b dkx * Ar^{k) + ® dky • Ay{k) (11) 

Jc Jc 

Since the energy contour under the influence of warp¬ 
ing is no longer a circle but has a dependence on (k,9), 




(16a) 

(16b) 


























3 


The Hall conductivity2ii2^ which is a^y = 777 ^ {jHaii,x) 

Ohiy 

(E X B) = 


and evaluating the curl such that 
gives 


dE, 


^ xy — 





(16c) 


Eq. IlScI by the factor— (1 — cosO) and integrating over 
all angles 


1 

r 


p r T{k,k'){i 

Jo 


cosO) d9 


ttP (?iw/fc)^+4A^ 
^ + (hvfk)'^'^ 


(19) 


where f{k) is the electron distribution of the given band. 
The Berry supported anomalous velocity^^^ in Eq. I16cl 
gives rise to a orbital magnetization(H fic) dependent 
Hall effect which flows transverse to the electric field. 
The anomalous velocity ceases to exist for a zero-gap 
system as can be easily derived by setting A to zero in 
Eq. [TU] and evaluating the curl of the Berry vector po¬ 
tential. In presence of a proximity induced band gap, 
further augmented by warping, the anomalous velocity is 
finite, a consequence of which is the intrinsic anomalous 
Hall conductivity. This anomalous Hall conductivity is 
therefore band gap and warping dependent. 

B. Momentum relaxation time 
The reciprocal of the relaxation time by solving the 
Boltzmann equation 

1 f d^k' 

-=I^J S,) \x,,' I (1 - COS 0 ) (17) 


where Xkk' — l(^ l^)PC j. The additional spin¬ 
scattering factor— ( ^s,s J takes in to account the helical 
spin structure of the TI surface states. If the external 
magnetic field is sufficiently large, such that the spins are 
aligned parallel to it, ( (^s,s 'j can be set to unity. /3 is a 
constant determined from nature of the scattering source. 
The wave vector scattering matrix can be computed as 
follows 


T{k,k) = {X+{k)exp{ie) A_(fc)) P (18a) 

where A± is given by Eq. I5bl Expanding the inner prod¬ 
uct gives 


T{k,k) 


1 

4 


[X\cos9 + A?_)^ -I- X\sin9 


Simplifying Eq. I18b[ one obtains 


(18b) 


Setting A to zero for a pristine 3D-topological insulator, 

the relaxation time reduces to — = —. The ratio of the 

T 4 


relaxation time is therefore 


Tt (fiu/fc)^-b 4A^ 
Tnt (hvfk)'^ + 

1 -b 4 k 
1 -b K 


( 20 ) 


where (Tnt) denotes the individual relaxation time for 

A2 

a topological (trivial) insulator and k = - If k ^ 

(hvfkf 

1 when a large magnetic field is impressed, the ratio is 
approximately equal to four. The scattering time for a 
topological insulator with a band gap is therefore reduced 
by a factor of four over its zero-gap counterpart. 

When the warping term is explicitly included, the band 
gap A must be modified to IX^arp = A -b tPXcosWk^. 
Using Aujarp, the scattering time expression in Eq. [TOl 
changes to 


1 2 (^A + XhPcos39k^)^ (1 — cos9) + {hvfk)^ sivPd 

Jo 2 ^{A + XhPcos39k^)^ + {hvfk)'^'^ 

( 21 ) 

The integral in Eq.[2T]is numerically evaluated in Section 

HI. 

The scattering time expression derived in Eq. can 
be inserted in a Boltzmann equation and solved under 
the relaxation time approximation to determine the re¬ 
sponse of the two-dimensional electronic system through 
the magneto-conductivity tensors axx and axy The 
magneto-conductivity tensors are evaluated with an ex¬ 
ternal magnetic field directed along the 2 :-axis and an 
electric field that is confined to the x-y plane. The elec¬ 
tric current in terms of Boltzmann transport equation 
per spin is therefore 


T{k,k) 


2A^ -b {hvfk)^ (1 -b cos9) 

2(A2-b(%fc)^) 


(18c) 


In deriving Eq. I18cl it is assumed that scattering takes 
place between two equi-energetic states such that |A:| = 
\k \. Further, the states are assumed to be of positive 
helicity (conduction band for a topological insulator) and 
exp{i9) is set to unity for the initial wave function(0 = 0). 
The relaxation time can now be evaluated by multiplying 


J = ^ / d^kvkSfk (22) 

where Sfk is the deviation from the Fermi-Dirac distri¬ 
bution, Vk = Vf {cos9,sin9) and 9 = tan~^ 

the Eermz-velocity. The Boltzmann equation assuming 
no spatial variation is 


df{r,k,t) M _ df{r,k,t) 

dt dt dk 


(23a) 
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where f{r,k^t) is the distribution function written as a 
sum of the equilibrium distribution and deviation under 
an electric and magnetic field f(r,k,t) = fo(r,k,t) + 

Sf{r, k, t). Since — = —g (E + v x B), Eq. I23al can be 
at 

written as 


J2-qf■E.Vk^ + i^rxB) 


dk 


= 0 (23b) 


dffr,k,t) 6f 

In writing Eq. I23b[ - —— is approximated as- 

ot T 

using the relaxation time approximation. E is the electric 
field on surface of the TI and r is carrier relaxation time. 
Solving Eq. I23bl and inserting Sf in Eq. [HI yields the 
magneto-conductivity tensors 


TtH"^ 1 -I- 


(24a) 


_ e^k/l 

TTh^ 1 + OJ^T^ 


(24b) 



FIG. 1: Dispersion of a 20.0 nm thick Bi 2 Se 3 topological in¬ 
sulator slab around the Dirac point (Fig. [1^) at 0.02 eV. The 
dispersion has a finite band gap (Fig. (TJd) when a ferromag- 
net is coated on the surfaces. The ferromagnet is assumed to 
possess an out of plane(«-directed) magnetization axis. The 
magnetization strength is an exchange energy set to 30.0 meV 


setting A to zero in Eq. 1141 Further, if the warping con¬ 
tribution is chosen such that it is equal in magnitude to 
the surface state energy because of the linear Hamilto¬ 
nian, the Berry phase expression changes to 


The ratio of the longitudinal and transverse conductivity 
tensors for a pristine and gapped topological insulator is 
therefore 


^ ^ Tnt (1 + 


and 


^ Tut (1 + 

Tf{l+UJ^T^t) 


(25a) 


(25b) 


where Wc = ev^Bje is the cyclotron frequency for Dirac 
fermions at a certain energy e. 

III. Results 

The presence of a non-topological component in a 
gapped 3D-TI (Fig. |T|) is numerically evaluated using ex¬ 
pressions derived in Sec II. The surface state of the TI 
is assumed to possess a finite band gap brought about 
most commonly through the proximity effect by interfac¬ 
ing with a ferromagnet that has an out-of-plane magne¬ 
tization component. 

A. Numerical evaluation of Berry phase with 
warping 

The combined influence of warping and a finite band 
gap on the Berry phase is computed numerically us¬ 
ing Eq. [131 The warping term /i^A(meVA^) and the 
band gap A are varied while |fc| is held constant at 0.03 
1/A. The accumulated Berry phase for a surface elec¬ 
tron with positive spin helicity in presence of proximity 
induced magnetic field and higher order warping terms 
is shown in Fig. [5J As the strength of the warping term 
and band gap increases, the deviation from the topolog¬ 
ically determined value of tt is more pronounced. The 
non-topological contribution arising purely on account of 
warping in absence of a band gap can be computed by 


which integrates to 1.2927r for a positive helicity elec¬ 
tron. This condition can be fulfilled when hvfk = 
Xh^k^cos{39). Setting cos (39) to 0.5, the warping term32 
to 200 evA^, and Vf = 5 x 10® m/s, gives k = 0.18 A. 
This roughly corresponds to an energy equal to 0.4 ev 
for the surface states. This in a way also determines a 
cut-off energy beyond which the warped Hamiltonian is 
the dominant component and outweighs the linear con¬ 
tribution. 

A significant manifestation of the increase in Berry 
phase can be seen as a way to enhance the Berry cur¬ 
vature (Eq. [T51) which is the geometric analog of a real 
magnetic field. Figure [5] shows the dependence of the di¬ 
mensionless ratio of the Berry curvature to square of the 

magnetic length ^ function of an exter¬ 

nal magnetic field. The magnetic length is approximated 
as 25.0 nm/'x/S. The band gap split is calculated using 
the Zeeman splitting where the ^-factor is set 

to 20 for electrons on surface of a TI. The Eermi-velocity 
for surface electrons in Bi 2 Te 3 is 5 x 10® m/s and the 
fc-vector is 0.3 1/A. The Berry curvature in Eg. [T5l does 
not include the warping term, which has been ignored to 
write a more compact expression. 

B. Momentum relaxation time 
The ratio of momentum relaxation times for the case 
of an ungapped topological insulator to a gapped sam- 

A2 . ^ , 

pie is roughly four times if k = - ^ in Eq. [20] is 

(hvfk) 

significantly larger than unity. The gap (A = gfisB) is 
computed by assuming an externally impressed magnetic 
field. Figure 0] shows the relaxation time ratio Tt/Tnt as 
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FIG. 2: The Berry phase ( 7 ) magnitude plotted as function of 
the warping strength. 7 increases with greater warping thus 
magnifying the non-topological contribution to the standard 
value of TT. The overall phase is shown for two values of the 
band gap indicated on the sub-plots. As warping increases, 
the band gap contribution diminishes and the curve flattens 
out. 



Magnetic field (T) 


FIG. 3: The ratio of Berry phase induced orbital magnetic 
field (Berry curvature) to magnetic length is plotted against 
an externally applied magnetic field. The Berry curvature is 
augmented as the magnetic field strength increases. Warping 
effects are ignored in this calculation. 

a function of the warping strength. The corresponding 
ratio for a magnetic field equal to 50 T without warping 
is also indicated as a constant. For relatively small values 
of the magnetic field, the linearly dispersing surface-state 
energy hvfk dominates the Zeeman splitting term. For 
instance, with B = 50 T along the 2 -axis, and a ^-factor 
equal to 20, the Zeeman-induced band gap is roughly 



Warping Constant (ev A^) 


FIG. 4: The relaxation time is computed for various warping 
strengths. At zero warping strength, the relaxation times for a 
finite band gap topological insulator coincide. The relaxation 
time is significantly altered due to a higher warping of the 
bands. 


0.0579 eV while the surface energy contribution at k- 
point = 0.3 1/A is 6.2 eV. The surface energy therefore 
masks the Zeeman split and brings the ratio of scattering 
times for a topological and trivial insulator close to unity. 
For /^-points sufficiently close to the Dirac point, which 
is a low energy case, for example, at fc = 5e-3 1/A, the 
ratio changes to 1.71. The overall validity of a surface 
state energy contribution limited to a linear Hamiltonian 
holds only when the fc-point is chosen close to the charge 
neutral Dirac point at F. The ratio of scattering times for 
points in momentum-space farther away from F is more 
correctly represented through a numerical evaluation of 
Eq. 1211 The warping term usually dominates the band 
gap split induced by the magnetic field and a marked 
deviation in the ratio of scattering times from the non- 
warped case is seen in Fig. 2] A higher warping strength 
significantly shifts the ratio curve away from values com¬ 
parable to unity for cases where the topological insulator 
is pristine or the higher-order k terms have limited con¬ 
tribution. 

The ratio of longitudinal and transverse conductivity 
tensors without warping is evaluated by using Eq. I25al 
and Eq. I25bl under a sufficiently strong magnetic field of 
10 T at 0.2 eV. The cyclotron frequency is 1.25 GHz; 
the ratio cr”*/cr”‘ by approximating 1 -f as 

is Tt/Tnt- The transverse conductivity under these con¬ 
ditions is almost close to unity. As mentioned above, 
at points in proximity to F, such as fc = 5e-3 1/A, the 
longitudinal conductivity tensor ratio is 1.71. The lon¬ 
gitudinal conductivity tensor therefore exhibits the same 
behaviour as the scattering times noted above. 

IV. Conclusion 

The influence of time reversal symmetry(TRS) break¬ 
ing magnetic field that opens a band gap at the surface 
and higher order terms defined by a warping of the eigen 









6 


energy spectrum substantially alter the topological Berry 
phase of (2n + 1) tt and any phenomenon that depends 
on the relaxation rate of the surface carriers. The Berry 
phase from the topologically determined value of tt gets 
an additional component that depends on the band gap 
split (TRS breaking) and the warping strength. For a 
given band gap, the Berry phase increases with the warp¬ 
ing strength. The altered Berry phase for a finite gap 
topological insulator described by a warped Hamiltonian 
also gives rise to anomalous Hall velocity through orbital 
magentization. The relaxation time within the Boltz¬ 
mann formalism (assuming spins aligned to the external 
magnetic field) is also warping dependent and at higher 
values offsets the magnetic field splitting. More realis¬ 
tic cases of time relaxation expression can be derived for 


charged impurities with screening and for the Kondo- 
effec t^°i^^ which involves spin-flip scattering by magnetic 
impurities. The spin scattering factor which has been 
chosen as unity to simplify calculations must be set as a 
scattering angle dependent quantity in a future work. 
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